Hyperbolicity of the Kidder-Scheel-Teukolsky formulation of Einstein's equations 
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We show that the Kidder-Scheel-Teukolsky family of hyperbolic formulations of the 3+1 evolution 
equations of general relativity remains hyperbolic when coupled to a recently proposed modified 
version of the Bona-Masso slicing condition. 
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I. INTRODUCTION 

The Cauchy problem for general relativity has re- 
ceived renewed interest in the last few years. To a 
large extent, this interest has been motivated by the 
realization that the mathematical structure of the evo- 
lution equations can have a direct impact on the sta- 
bility of numerical simulations. Research has concen- 
trated in developing strongly, or even symmetric hyper- 
bolic formulations of the evolution equations of general 
relativity, as such systems can be shown to be well- 
posed [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. 
This well-posedness implies that one can find numeri- 
cal discretizations that are stable in the sense that the 
growth of errors is bounded [16]. 

A related problem to that of finding well-posed systems 
of evolution equations is the problem of finding well be- 
haved coordinate systems. In a 3+1 formulation, this 
problem reduces to choosing conditions that determine 
the so-called "gauge" quantities, that is, the lapse func- 
tion and shift vector. The lapse function determines the 
slicing of the 4-dimcnsional spacetime into 3-dimensional 
spatial hypersurfaces, and the shift vector relates the 
spatial coordinate systems of nearby hypersurfaces. Our 
group has recently concentrated in studying slicing con- 
ditions that can be written as hyperbolic equations for a 
time function T whose level surfaces correspond to the 
members of the foliation [17, 18]. In Ref. [17] we concen- 
trated in the so-called Bona-Masso (BM) family of slicing 
conditions [3] and studied under which circumstances it 
avoids different types of pathological behaviors, while in 
Ref. [18] we proposed a modified version of the BM slicing 
condition that is well adapted to the evolution of static 
or stationary spacetimes and to the use of a densitized 
lapse as the fundamental variable. 

Whenever one proposes a new gauge condition, the is- 
sue arises of studying if such a condition affects the well- 
posedness of the system of evolution equations as a whole. 



Such an analysis, for example, has been carried out by 
Sarbach and Tiglio [19] for a generalization of the BM 
condition and more recently by Lindblom and Scheel [20] 
for another generalization of the BM condition coupled to 
a 'T-driver" shift condition [21]. In both these cases the 
analysis was done using multi-parameter first order for- 
mulations of the Einstein evolution equations. Here we 
will consider the Kidder-Scheel-Teukolsky (KST) formu- 
lation [15] coupled to the modified BM slicing condition 
studied in [18]. 

This paper is organized as follows. In Sec. II we in- 
troduce briefly the BM slicing condition and its modi- 
fied form. Section III describes the KST formulation of 
the Einstein evolution equations. In Sec. IV we analyze 
the hyperbolicity of the coupled system of KST evolu- 
tion equations plus modified BM slicing condition. We 
conclude in Sec. V. 



II. THE MODIFIED BONA-MASSO SLICING 
CONDITION 

The BM family of slicing conditions [3] is well known 
and has been discussed extensively in the literature (see 
for example [17, 22] and references therein). This slic- 
ing condition asks for the lapse function to satisfy the 
following evolution equation 



-a = (d t -C (i )a = -a 2 f(a)K, 



(2.1) 
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with £p the Lie derivative with respect to the shift vec- 
tor (3 l , K the trace of the extrinsic curvature and f(a) a 
positive but otherwise arbitrary function of a. This con- 
dition can be shown to be hyperbolic in the sense that it 
is equivalent to asking for the time function T to satisfy 
a generalized wave equation. 

In a recent paper [18], we have proposed a modified 
version of condition (2.1) that keeps many of its impor- 
tant properties but is at the same time well adapted to 
the evolution of static or stationary spacetimes and also 
to the use of a densitized lapse as a fundamental vari- 
able. We believe that having a slicing condition that is 
compatible with a static solution is a necessary require- 
ment if one looks for symmetry seeking coordinates of 
the type discussed by Gundlach and Garfinkle [23] and 
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by Brady et.al [24], that will be able to find the Killing 
fields that static (or stationary) spacetimes have, or the 
approximate Killing fields that many interesting astro- 
physical systems will have at late times. This modified 
BM slicing condition has the form 



d t a = -af(a) (aK - V,/?) , 



(2.2) 



with Vj is the 3-covariant derivative associated with f/^ . 
One can show that this condition can also be obtained 
from a generalized wave equation for the time function 
T and is hence also hyperbolic independently of the Ein- 
stein equations. 



III. THE KST FAMILY OF FORMULATIONS OF 
THE EINSTEIN EVOLUTION EQUATIONS 

The KST family of formulations of the Einstein evo- 
lution equations is a multi-parameter, fully first order, 
system of equations for 30 independent dynamical vari- 
ables {gij, Kijjdf-ij}, where is the spatial metric, 
the extrinsic curvature, and dkij '■= dkgij- Notice that 
the definition of the dkij is used only for obtaining initial 
data, the d's are then promoted to independent variables 
and their definition in terms of derivatives of the g's then 
becomes a constraint. 

If we define do = (d t — Cp)/a, the system of evolution 
equations in vacuum can be written as: 



= -2^-, 

= Rij -(V i V j a)/a-2K im K r J 
+ KK l3 + 7 9ij C + C g ah C a(mh , 



d d, 



kij 



-2 d k - 2K lj d k In a 
+ il 9k(iCj) + X 9ijCk , 
where {7, £, 77, \} are f ree parameters and 



C 

Clkij 



(R - K ab K ab + K 2 ) /2 . 
W a K m - WiK, 

d[ldf~]ij , 



(3.1) 
(3.2) 
(3.3) 

(3.4) 
(3.5) 
(3.6) 
(3.7) 



are constraints of the system (the first two are the Hamil- 
tonian and momentum constraints, and the last two are 
consistency constraints). Notice that since the dkij are 
not components of a tensor, their Lie derivative with re- 
spect to (3 l should be understood as 

Cpdkij = [3 a d a dki] + d aij dk(3 a 

+ 2d ka(l d j) p a + 2g a{l d ]) dkp a . (3.8) 

The Ricci tensor that appears in the evolution 
equation for is written in terms of the d's as 
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dad 



(lb 



(—d a dbij + d a d(ij) b + d(id\ ab \j) 



1 

+ 2 



di ab d jab + (d k - 2b k ) r 



with dk ■= g v dkij,bk := g tJ dijk and Tj k the Christoffel 
symbols associated with g^. It is important to mention 
that the system of equations above is not the most gen- 
eral form of the KST system which has 12 free param- 
eters. Here we have considered only the 4 parameters 
that are related to constraint terms and ignored the 7 
parameters that redefine the independent variables and 
the parameter related to the weight of the prescribed den- 
sitized lapse which we will replace with our modified BM 
slicing condition. 

In the original analysis of KST, the system of equa- 
tions (3.1)-(3.3) was shown to be strongly or even sym- 
metric hyperbolic for certain regions of the parameter 
space {jXtVjXIi with the lapse replaced by a "densi- 
tized lapse" q given by 



q := Hg-°a) 



(3.10) 



with g the determinant of g^ , and a positive (with a pre- 
ferred value of 1/2). The densitized lapse q was assumed 
to be a prescribed, i.e. a priori known, function of space 
and time. This condition was later relaxed by Sarbach 
and Tiglio in [19] where the lapse was instead taken to 
be an arbitrary function of g such that 



Ceff := — d q a > . 
a 



(3.11) 



IV. HYPERBOLICITY OF THE KST 
FORMULATION COUPLED TO THE MODIFIED 
BM CONDITION 

We start from the modified BM slicing condition (2.2) 
which we rewrite as 



d t a = -af{a)T 



with 



T:=aK- V m /3 m . 
We now define the first order quantity: 

di In a 

From Eq. (4.1) one can easily show that 
d t A t = -3{T . 



(4.1) 
(4.2) 

(4.3) 

(4.4) 



On the other hand, the derivatives of a that appear 
in the evolution equation for given in the previous 
section, Eq. (3.2), can be written in terms of Ai as 

= f [d {t A j} + (/ + a/') AiAj T%A k ] , (4.5) 

where we have used the fact that diAj is symmetric. 
Notice now that from the evolution equation for g^, 
Eq. (3.1), one can also find that 



1 im 1 jk 7 



(3.9) 



d t g = -2gT, 



(4.6) 
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which implies that 



dtD, = -2 d(T 



(4.7) 



with Di := dilng. Comparing equations (4.4) and (4.7) 
we find 



d t A t = - dtD, . 



(4.8) 



Now, from the definition of dkij, we should have 
Di = di, with di as defined in the previous section. How- 
ever, since in the KST formulation the evolution equa- 
tions for the dkij are modified by adding multiples of 
constraints to them, we will generally have d t Di ^ d t di. 
Because of this, we propose to modify the evolution equa- 
tion for Ai in the following way 

dtAi = —diT + Fi(C,Ck,Ckim,Ckimn) ■ (4.9) 

From the evolution equation (3.3) for dkij, one can find 
after some algebra 



d t d t = -2 d{T + a(r] + 3 X )C, 
+ 2aK ah C tab + C ma a d^ 
+ f3 m d m Ci a a , 



(4.10) 



which means that if we take 



2Fi = a( V + 3 x )Ci + 2aK ab C lab 

+ C ma a d t p m + (3 m d m C ta a , (4.11) 



then we will always have 



1 



d t Ai = - d t di . 



(4.12) 

The last equation allows us to define the quantities 

Q t := Ai - d t /2 . (4.13) 

These quantities are then such that d t Qi — 0, that is, 
they are non-dynamical. 

Another way to introduce the Qi is the following: From 
the modified BM condition and the evolution equation for 
gij it is easy to show that 



d t a _ dtg_ 
af ~ 2g ' 



(4.14) 



which one can easily integrate to find 

gV^Htf) ex P /^, (4.15) 

with H(x l ) an arbitrary time-independent function. This 
shows that if we define 



q : 



:= In a- 1 ' 2 



g 1 exp 



da \ 



(4.16) 



then we will have dtq = 0. Notice that the q defined 
above is just the generalization of the densitized lapse 



defined in (3.10) for the case / ^ 1. One can now show 
that the Qi defined through (4.13) are precisely such that 
Qi = diq, and since q is time independent, then so are 
the Qi. 

Having introduced the non-dynamical quantities Qi, 
we can rewrite the derivatives of Ai appearing in the 
evolution equation of Kij through the term (4.5) in terms 
of derivatives of Qi and di. Since the Qi do not evolve, 
they can be considered as source terms. In this way, the 
system of evolution equations for and dkij becomes 

d Kij ~ \g ab [-d a d blJ + (l + ()d a d {ll)b 

+ (1 - C) d {i d\ ab \j) - (1 + /) d {l dj )ab 

+ 7 9 l] g kl da(dkib-d bk i)} , (4.17) 

dodkij ~ -2d k Kij +rig Hl g ah {d\ a \K j)b - d 3) K ab ) 

+ X 9^9 ab (d a K kb - d k K ab ) , (4.18) 

where the symbol ~ means equal up to principal part. 
The system above is exactly the same as the one pre- 
sented by Sarbach and Tiglio in [19] with the replacement 
(T c ff = J/2. The hyperbolicity analysis of that reference 
then follows directly. In particular, the non-zero eigen- 
values of the system become 



Ai = /, 



(4.19) 



A 2 = 1 + X- -(l + C>? + 7(2-»? + 2x), (4.20) 
As = |x+|(l-C)»7- J(l + /)fa + 3x), (4.21) 



A 4 



1 



(4.22) 



There are 12 eigenvectors associated with these non-zero 
eigenvalues: two with both Ai and A 2 , and four with 
both A3 and A4. There are 12 more eigenvectors with 
eigenvalue zero. The system can be shown to be strongly 
hyperbolic if 

Xj > 0, for j = 1,2,3, 
A 3 = |(3Ai + 1) if Ai =A 2 . 

The associated characteristic speeds are given simply by 
vf = i(Ai) 1 / 2 . In particular, we obtain vf = i/ 1 / 2 , 
which agrees with the expected result for the BM slicing 
condition. 



V. DISCUSSION 

We have studied the hyperbolicity of the KST family 
of formulations of the Einstein evolution equations cou- 
pled to a recently proposed modified BM slicing condi- 
tion. We find that the modified BM condition allows one 
to construct a non-dynamical function q that generalizes 
the densitized lapse to the case when the function f(a) 
defining the slicing is different from 1. From this non- 
dynamical quantity one can construct three first order 
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non-evolving quantities Qi := diq that can be used to re- 
place the spatial derivatives of the lapse in the evolution 
equation of the extrinsic curvature Kij. By doing this 
we are able to reduce the system of evolution equations 
to one previously analyzed by Sarbach and Tiglio, which 
allows us to show that the coupled KST plus modified 
BM slicing condition system remains strongly hyperbolic 
in the same circumstances as before, and also to identify 
directly the characteristic speeds. 
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